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ABSTRACT 

Kilogauss-strength magnetic fields are often observed in intergranular lanes at the 
photosphere in the quiet Sun. Such fields are stronger than the equipartition field B e , 
corresponding to a magnetic energy density that matches the kinetic energy density 
of photospheric convection, and comparable with the field B p that exerts a magnetic 
pressure equal to the ambient gas pressure. We present an idealised numerical model 
of three-dimensional compressible magnetoconvection at the photosphere, for a range 
of values of the magnetic Reynolds number. In the absence of a magnetic field, the 
convection is highly supercritical and is characterised by a pattern of vigorous, time- 
dependent, "granular" motions. When a weak magnetic field is imposed upon the con- 
vection, magnetic flux is swept into the convective downflows where it forms localised 
concentrations. Unless this process is significantly inhibited by magnetic diffusion, the 
resulting fields are often much greater than B e , and the high magnetic pressure in 
these flux elements leads to their being partially evacuated. Some of these flux el- 
ements contain ultra-intense magnetic fields that are significantly greater than B p . 
Such fields are contained by a combination of the thermal pressure of the gas and the 
dynamic pressure of the convective motion, and they are constantly evolving. These 
ultra-intense fields develop owing to nonlinear interactions between magnetic fields 
and convection; they cannot be explained in terms of "convective collapse" within a 
thin flux tube that remains in overall pressure equilibrium with its surroundings. 

Key words: convection - (magnetohydrodynamics:) MHD - Sun: granulation - Sun: 
magnetic fields - Sun: photosphere 



1 INTRODUCTION 

High resolution observations of the solar photosphere have 
greatly enhanced our understanding of the complex interac- 
tions between magnetic fields and turbulent convective mo- 
tions. Modern ground-based telescopes, such as the SST on 
La Palma, and space-borne instruments, such as the Solar 
Optical Telescope (SOT) on the recently-launched Hinode 
satellite, can now resolve the fine structure of photospheric 
magnetic fields. Such high resolution observations have al- 
ready provided new insights into the magnetic field struc- 
tures that can be found in sunspot penumbrae (Scharmer et 
al. 2002) , plage regions (Berger et al. 2004) and the quiet Sun 
(see, for example, Centeno et al. 2007; Rezaei et al. 2007). 
In the quiet Sun, localised concentrations of intense verti- 
cal magnetic flux often accumulate in the convective down- 
flows (Lin & Rimmele 1999; Domfnguez Cerdefia, Kneer & 
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Sanchez Almeida 2003; Centeno et al. 2007). These localised 
flux concentrations typically show up as bright points in G- 
band images of the photospheric granulation (see, for ex- 
ample, Berger & Title 1996, 2001). Direct measurements 
indicate that the peak field strength within these small- 
scale magnetic features is often well in excess of a kilogauss 
(see, for example, Grossmann-Doerth, Keller & Schiissler 
1996; Dommguez Cerdefia et al. 2003). Although kilogauss- 
strength magnetic features occupy only a small fraction of 
the quiet solar surface, these localised flux concentrations 
contain a significant proportion of the total (unsigned) quiet 
Sun magnetic flux (Dommguez Cerdefia, Sanchez Almeida 
& Kneer 2006; Sanchez Almeida 2007). 

The origin of quiet Sun magnetic fields remains an open 
question. Either these magnetic fields are generated locally 
at the photosphere, as the result of small-scale convectively- 
driven dynamo action (Cattaneo 1999; Cattaneo & Hughes 
2006; Vogler & Schiissler 2007), or they were originally gen- 
erated elsewhere in the solar interior and are simply being 
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re-processed and amplified by the local convective motions. 
What is clear, however, is that the granular convective up- 
wellings at the solar photosphere will tend to expel magnetic 
flux, causing it to accumulate in the convective downflows 
(Proctor & Weiss 1982). This explains the observed associa- 
tion between quiet Sun magnetic fields and the intergranular 
lanes. However, the peak field strength that is measured in 
these regions is more difficult to explain. Estimates of the ki- 
netic energy density of the granular convection in the solar 
photosphere suggest that magnetic field strengths greater 
than approximately 400G would exceed the "equipartition" 
value, B e , at which the magnetic energy density of these 
localised fields is comparable to the granular convective ki- 
netic energy density (see, for example, Galloway, Proctor 
& Weiss 1977). Since the magnetic energy density is pro- 
portional to the square of the magnetic field strength, the 
observed kilogauss-strength fields in the quiet Sun have a 
magnetic energy density that is an order of magnitude larger 
than equipartition. In fact, a better estimate for the observed 
field strengths is given by B p , which corresponds to the field 
strength at which the local magnetic pressure balances the 
ambient granular gas pressure. 

In theoretical studies of Boussinesq magnetoconvection 
(see e.g. Galloway et al. 1977; Galloway, Proctor & Weiss 
1978; Galloway & Moore 1979), the process of magnetic flux 
expulsion leads naturally to the formation of localised flux 
concentrations at the edges of convective cells. For axisym- 
metric converging flows, the kinematic flux concentration of 
weak magnetic fields produces peak fields that scale linearly 
with the magnetic Reynolds number, Rm, whilst the width 
of the resulting magnetic feature is proportional to RmT 1 ^ 2 . 
However, for stronger fields this amplification process is lim- 
ited by dynamical effects. Since the magnetic pressure does 
not play a role in Boussinesq magnetoconvection, any mag- 
netic feedback is due to the non-conservative component of 
the Lorentz force, which results from the magnetic curva- 
ture. In the dynamical regime, the peak field can only exceed 
B e for large values of the magnetic Prandtl number (the ra- 
tio of the kinematic viscosity of the fluid to its magnetic dif- 
fusivity) . Formally, this scaling can only be verified for small 
Reynolds numbers, Re, although it appears not to be re- 
stricted to this parameter regime (Cattaneo 1999). Kerswell 
& Childress (1992) also found a similar magnetic Prandtl 
number dependence in their idealised boundary-layer study 
of the equilibrium of a thin flux tube in steady compress- 
ible convection (see also Cameron & Galloway 2005). Un- 
fortunately, estimates of the magnetic Prandtl number in 
the solar photosphere indicate that this diflusivity ratio is 
extremely small (see, for example, Ossendrijver 2003). So 
although these comparatively idealised models can explain 
the formation of highly localised magnetic features at large 
Rm, they are unable to account for the appearance of super- 
equipartition, kilogauss-strength magnetic fields in the quiet 
Sun. 

Since the peak magnetic field in kilogauss-strength mag- 
netic features is comparable to B p , the magnetic pressure 
is bound to play a significant role in the local dynamics. 
If a magnetic feature is assumed to be in pressure balance 
with its non-magnetic surroundings, the internal gas pres- 
sure must be smaller than that of the surrounding fluid in 
order to compensate for the large magnetic pressure. Un- 
less these features are much cooler than their surroundings, 



this reduction in gas pressure is most easily achieved if these 
magnetic regions are at least partially evacuated (see, for ex- 
ample, Proctor 1983; Proctor & Weiss 1984). This effect is 
absent from Boussinesq models, which neglect the effects of 
compressibility, so it is inevitable that such models underes- 
timate the peak fields that can be generated by the process of 
flux concentration at the solar photosphere. When seeking to 
explain the formation of kilogauss-strength magnetic fields 
in the quiet Sun, it is therefore important to consider the 
effects of compressibility (see e.g. Hughes & Proctor 1988). 

The most widely studied compressible model of mag- 
netic field amplification in the quiet Sun is the "convective 
collapse" instability of thin magnetic flux tubes (Webb & 
Roberts 1978; Spruit 1979; Spruit & Zweibel 1979; Unno 
& Ando 1979). Convective collapse models consider the lin- 
ear stability of a thin vertical magnetic flux tube embed- 
ded in a non-magnetic, superadiabatically stratified atmo- 
sphere. There are initially no convective motions along the 
tube, which is assumed to be in pressure balance with its 
surroundings. Provided that the initial magnetic field is not 
too strong, the flux tube is subject to a convective insta- 
bility that drains plasma vertically out of the tube. As the 
plasma drains out of the tube, the local pressure balance 
implies that the tube must "collapse" to form a narrower 
(and therefore more intense) magnetic flux concentration. 
This instability will continue to operate until the magnetic 
field is strong enough to suppress convective motions. This 
model certainly represents a plausible mechanism for the 
formation of kilogauss-strength magnetic fields (the only re- 
striction being that the peak field, B max ^ B p ). Having 
said that, this model has its limitations (Hughes & Proc- 
tor 1988; Thomas & Weiss 1992, 2008). In particular, the 
thin flux tube approximation does not appear to be consis- 
tent with photospheric observations (see, e.g. Berger & Title 
1996; Berger et al. 2004): even in the quiet Sun, the observed 
magnetic regions seem to form a constantly-evolving "mag- 
netic fluid" (see, for example, Thomas & Weiss 2008), and 
so cannot simply be regarded as a collection of discrete col- 
lapsed flux tubes. In addition, the static equilibrium for the 
instability is rarely achievable in photospheric magnetocon- 
vection. Finally, by reducing the horizontal component of 
the momentum equation to a simple pressure balance, the 
model not only ignores the possible dynamical influence of 
the surrounding convective motions, but also neglects the 
possible effects of magnetic curvature, which may play an 
important role in ensuring the equilibrium of more exten- 
sive magnetic features (see, for example, Vogler et al. 2005, 
where curvature effects help to support a peak field, B ma x, 
in excess of -Bp). 

More realistic models of this process can only be studied 
by carrying out large-scale numerical simulations of convec- 
tive magnetic flux intensification. In order to make detailed 
comparisons between numerical simulations and spectral ob- 
servations, it is necessary to take account of effects such as 
partial ionisation and radiative transfer (see, e.g., Nordlund 
1982; Vogler et al. 2005). This approach has successfully sim- 
ulated magnetic fields and spectral features that are similar 
to those observed in plages (Grossmann-Doerth, Schussler & 
Steiner 1998; Keller et al. 2004; Carlsson et al. 2004; Vogler 
et al. 2005), in sunspot umbrae (Schussler & Vogler 2006), 
and in the quiet Sun (Khomenko et al. 2005; Stein & Nord- 
lund 2006). More idealised models of photospheric magne- 
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toconvection focus entirely upon the interactions between 
compressible convection and magnetic fields (e.g. Hurlburt 
& Toomre f988; Matthews, Proctor & Weiss 1995; Weiss 
et al. 1996; Rucklidge et al. 2000; Weiss, Proctor & Brown- 
john 2002; Bushby & Houghton 2005). This complementary 
approach lends itself more easily to systematic surveys of 
parameter space, and has had considerable success in qual- 
itatively reproducing solar-like behaviour. 

In this paper, we investigate the formation of super- 
equipartition, localised magnetic features in the quiet Sun. 
We carry out numerical simulations of an idealised model 
of photospheric magnetoconvection. Clearly it is not possi- 
ble to achieve solar-like values of the viscous and magnetic 
Reynolds numbers in these simulations (which are both very 
large at the solar photosphere, see e.g. Ossendrijver 2003). 
However the Reynolds numbers that are used here are large 
enough that these simulations (at least qualitatively) repro- 
duce most of the key physical features of flux intensifica- 
tion in photospheric magnetoconvection. The set-up of this 
model and the numerical results that are obtained from it 
are described in the next two Sections of the paper. In the 
final part, we discuss the relevance of these results to photo- 
spheric magnetic fields and compare our findings with pre- 
dictions of the simplified "convective collapse" models. 



main (which are denoted by To and po respectively). Defin- 
ing i?» to be the gas constant, the velocity, u, is scaled 
in terms of the isothermal sound speed at the top of the 
layer, (7?»To) 1//2 ; a natural scaling for time is therefore 
d/(i?*Tb) , which corresponds to an acoustic timescale. 
The parameter 8 denotes the (dimensionless) temperature 
difference between the upper and lower boundaries. In the 
absence of any motion the gas is a polytrope with a poly- 
tropic index m = gd/R t To9 — 1. Rather than relating the 
magnetic field strength to the Chandrasekhar number (see, 
for example, Weiss et al. 2002; Bushby & Houghton 2005), 
we also scale the Alfven speed at the top of the layer in 
terms of the sound speed. This implies that the appropri- 
ate non-dimensionalising scaling for the magnetic field, B, 
is given by (popoR^To) 1 ^ 2 . Having made these scalings, the 
governing equations for the density, the momentum density 
(pu), the magnetic field and the temperature are given by 

| = -V.(p U ), (1) 

^- t (pu) = -V(P+\B\ 2 /2)+e(m+l)pz (2) 
+V • (BB - puu + kot) , 



2 GOVERNING EQUATIONS AND 
NUMERICAL METHODS 

The model that is considered in this paper is an idealised, 
local representation of magnetoconvection in the quiet Sun. 
We solve the equations of three-dimensional compressible 
magnetohydrodynamics for a plane layer of electrically- 
conducting, perfect monatomic gas. The layer of gas is 
heated from below and cooled from above. We assume con- 
stant values for the gravitational acceleration g, the shear 
viscosity p, the magnetic diffusivity rj, the thermal conduc- 
tivity K, the magnetic permeability po and the specific heat 
capacities at constant density and pressure (c v and c p re- 
spectively). The axes of the chosen Cartesian frame of refer- 
ence are oriented so that the z-axis points vertically down- 
wards (parallel to the gravitational acceleration). Defining 
d to be the depth of the convective layer, the gas occupies 
the region ^ x,y ^ Ad and < z < d, which gives a wide 
Cartesian domain with a square horizontal cross-sectional 
area. Periodic boundary conditions are imposed in each of 
the horizontal directions. Idealised boundary conditions are 
imposed at z = (the upper surface) and z = d (the lower 
surface) . These bounding surfaces are held at fixed tempera- 
ture, and are assumed to be impermeable and stress-free. In 
addition, any magnetic field that is present is constrained to 
be vertical at the upper and lower boundaries. Similar mod- 
els have been considered in several previous studies (see, for 
example, Matthews et al. 1995; Rucklidge et al. 2000; Weiss 
et al. 2002; Bushby & Houghton 2005). 

It is convenient to formulate the governing equations for 
magnetoconvection in terms of non-dimensional variables. 
We adopt non-dimensionalising scalings that are similar to 
those described by Matthews et al. (1995). All lengths are 
scaled in terms of the depth of the Cartesian domain, d, 
whilst the temperature, T, and density, p, are both scaled 
in terms of their values at the upper surface of the do- 



— = V x (u x B - <oV x B) , V • B = 0, (3) 

^ = -w-VT-( 7 -1)TV-u+^Vt (4) 

+ (^/2 + ColVxBI 2 ). 

The components of the stress tensor, r are given by 

T% - du i | du i 2duk 5i ; ( 5 ) 
%3 dxj dxi 3 dx k l; " 

whilst the pressure, P is determined by the equation of state 
for a perfect gas 

P = P T. (6) 

These governing equations are characterised by sev- 
eral non-dimensional constants, including the Prandtl num- 
ber, a — pCp/K, the (non-dimensional) thermal diffusiv- 
ity, k = K/p dCp (TitTo) 1 / 2 , and the ratio of specific heats, 
7 = Cp/cv. The ratio of the magnetic diffusivity to the 
thermal diffusivity, f , is proportional to the fluid density p 
(and therefore increases with depth). At the top of the layer 
C = Co = rjpoCp/K. If all the other parameters are fixed, 
varying k is equivalent to varying the mid-layer Rayleigh 
number, 

Ra = (m + 1 - m 7 ) (1 + g/2) 2 — 1 (m ± . (7) 

This Rayleigh number measures the destabilising effects of 
a superadiabatic temperature gradient relative to the stabil- 
ising effects of diffusion. The parameters are all described in 
greater detail by Matthews et al. (1995). 

The simplest non-trivial equilibrium solution of these 
governing equations corresponds to a static polytrope with 
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a uniform magnetic field. Initially we restrict attention to the 
case in which B = 0. In this equilibrium solution, u = 0, 
T(z) = 1 + Oz and p(z) = (1 + 9z) m . Fixing 9 = 10 and 
m — 1 gives a highly stratified atmosphere in which the 
temperature and density both vary by an order of magni- 
tude across the layer. We consider a monatomic gas, there- 
fore 7 = 5/3. This implies that the m = 1 polytrope is 
superadiabatically stratified. For simplicity (and for ease of 
comparison with previous studies) we set the Prandtl num- 
ber equal to unity, i.e. a = 1. A range of values for is 
considered in this paper. Finally, the Rayleigh number is 
chosen to be Ra = 4.0 x 10 s . This value for Ra is more than 
two orders of magnitude larger than the critical value for 
the onset of convective instabilities. 

The model described above is investigated by carrying 
out large-scale numerical simulations. For this 4x4x 1 Carte- 
sian domain, we adopt a computational mesh of 256 x 256 x 
160 grid points. The horizontal extent of this computational 
domain is smaller than in some of our previous calculations 
(see, e.g. Bushby & Houghton 2005). This reduction in box 
size does influence the global structure of the convective mo- 
tions by eliminating mesoscale structures (Rucklidge et al. 
2000). However, the horizontal extent of the computational 
domain probably has little influence upon the localised pro- 
cess of flux concentration. The advantage of considering a 
smaller domain is that it allows us to carry out high reso- 
lution numerical simulations without requiring an excessive 
number of grid points. This enables us to model magne- 
toconvective behaviour more accurately at high Reynolds 
numbers, so this reduction in box size is a reasonable com- 
promise. We use a well tested code (see, e.g. Matthews et 
al. 1995) in which horizontal derivatives are evaluated in 
Fourier space, whilst vertical derivatives are evaluated us- 
ing fourth order finite differences. The time-stepping is car- 
ried out via an explicit 3rd order Adams-Bashforth scheme, 
with a variable time-step. The code is efficiently parallelised 
using MPI, and these simulations have made use of the 
Cambridge-Cranfield High Performance Computing Facility 
and the UKMHD Consortium machine based in St Andrews. 



3 RESULTS 

3.1 The initial state 

The starting point for these idealised simulations is fully 
developed non-magnetic convection. In order to generate 
such a convective state, we consider the equilibrium solu- 
tion corresponding to a static polytrope with no applied 
magnetic field, and introduce a small amplitude random 
temperature perturbation. This convectively-unstable con- 
figuration is allowed to evolve until it has relaxed to a sta- 
tistically steady hydrodynamic state, as illustrated in Fig. [1] 
which shows a snapshot of the temperature distribution in 
a horizontal plane near the upper surface of the compu- 
tational domain. In this granular pattern, bright regions 
correspond to broad warm upflows, whilst darker regions 
correspond to cool narrow downflows. The strength of the 
convection can be characterised by the mid-layer Reynolds 
number Re = PmidUrmsd/ fi, where U lms corresponds to the 
rms-velocity of the convection and p m id is the density at 
the mid-layer of the original static polytrope. The Reynolds 
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Figure 1. The temperature distribution for non-magnetic con- 
vection, in a horizontal plane near the upper surface of the Carte- 
sian domain, at (z/d) = 0.05. Temperature contours are evenly 
spaced between T = 1.4 and T = 2.6 (the same contour spacings 
are used for all temperature plots in this paper). Brighter regions 
correspond to warmer fluid, cooler areas of fluid arc represented 
by darker regions. 

number here is approximately Re = 150. As noted in the 
Introduction, realistic Reynolds numbers for photospheric 
convection are numerically unobtainable, but the Reynolds 
number is large enough that instructive results can be ob- 
tained in these idealised calculations. 

Having established this (purely hydrodynamic) convec- 
tive state, we then introduce a weak uniform vertical mag- 
netic field, Bqz, with Bo chosen so that the initial magnetic 
energy is approximately 0.1% of the kinetic energy. A weak 
imposed field of this form tends to favour the formation of 
highly localised magnetic features (see, for example, Bushby 
& Houghton 2005). In what follows, the time at which this 
magnetic field is introduced is denoted by t — 0. The subse- 
quent evolution of this magnetic field depends crucially upon 
the magnetic Reynolds number of the flow, Rm — Uruiad/r/. 
In non-dimensional terms, Rm oc ; thus it is possible 
to investigate a range of values of Rm simply by repeating 
the numerical experiment with different values of (o. The 
range investigated is 0.2 Co ^ 2.4 (120 > Rm > 10). 
Larger values of Rm correspond to less diffusive plasmas, 
and are therefore more relevant to photospheric magneto- 
convection, although (as with Re) it is not yet possible to 
carry out fully resolved simulations with realistic values of 
Rm for photospheric magnetoconvection. We therefore fo- 
cus initially upon the case of £o = 0.2 (Rm ~ 120). The 
effects of varying £o are discussed later in this Section. 

3.2 Results for Co = 0.2 

Since the initial magnetic field is comparatively weak, the 
Lorentz forces play a minor role in the very early stages of 
evolution in these numerical simulations. During this brief 
"kinematic" phase, diverging convective motions at the up- 
per surface of the computational domain rapidly expel mag- 
netic flux from the granular interiors. This process causes 
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Figure 2. The magnetic field distribution during the flux ex- 
pulsion phase, at t = 0.24. The shaded contours show the tem- 
perature distribution in a horizontal plane at (z/d) = 0.05 (as 
described in Fig. [TJ. The white contours denote constant values 
of the vertical component of the magnetic field. 



magnetic flux to accumulate in the convective downflows in 
the intergranular lanes. The strongest field concentrations 
tend to occur at the vertices between neighbouring gran- 
ules, due to the fact that any flows along the intergranu- 
lar lanes tend to converge upon these vertices. The process 
of flux expulsion and accumulation is illustrated in Fig. [2] 
which shows the distribution of the vertical component of 
the magnetic field during this flux expulsion phase. In this 
(relatively) high Rm regime, the width of the localised mag- 
netic features is comparable to the width of the narrow in- 
tergranular lanes. 

As these localised concentrations of magnetic flux form, 
the principle of flux conservation implies that the mag- 
netic energy density in these regions rapidly increases to 
the point where the Lorentz forces become dynamically sig- 
nificant. In the upper layers of the computational domain, 
most of the magnetic energy resides in the vertical com- 
ponent of the magnetic field, and the horizontal gradients 
in this field component are usually much larger than any 
variations in the vertical direction. This implies that the 
horizontal magnetic pressure gradient (i.e. the horizontal 
component of V [B 2 /2/^o]) is the dominant component of 
the Lorentz force. The effects of the magnetic pressure are 
most apparent near the surface, where the gas pressure is 
relatively small. The magnetic pressure gradient tends to 
inhibit the converging convective motions that are responsi- 
ble for driving the flux concentration process. However, the 
flux amplification process is not immediately suppressed, be- 
cause the total (i.e. gas plus magnetic) pressure increases 
more gradually than the magnetic pressure alone. This is 
due to the fact that the local convective downflows rapidly 
carry fluid away from the surface regions of the magnetic 
feature. These downflows are responsible for partially evac- 
uating the upper regions of the magnetic features, which in 
turn leads to a reduction in the local gas pressure. This re- 
duction in gas pressure (at least partially) compensates for 
the increased magnetic pressure. In the most intense mag- 



netic flux concentrations, the gas pressure drops to as little 
as a few percent of its initial value at the surface before the 
magnetic field becomes locally strong enough to suppress 
the vertical convective motions. A combination of the sup- 
pression of these vertical motions plus the effects of diffusion 
eventually halts this flux concentration process. 

To describe this process in a more quantitative fashion, 
it is useful not only to consider the variations of the gas 
pressure (-P ga s) and the magnetic pressure (P mag ), but also 
the dynamic pressure (Pdyn), which represents the dynami- 
cal influence of the convective motions. Following Hurlburt 
& Toomre (1988), we define -Pdyn = p| u | 2 (although note 
that other definitions have been used in some other stud- 
ies, e.g., Weiss et al. 1996). Whilst this expression does not 
correspond directly to a pressure term in the momentum 
equation (2), it does usefully quantify the vigour of the con- 
vective motions. Since the peak flow speeds at the surface 
are comparable to the local sound speed (any shocks be- 
ing smoothed out by the viscosity of the fluid), we would 
expect -Pdyn to play a significant dynamical role. Fig. 3 il- 
lustrates the time-evolution of these pressure distributions 
at the surface (where the most significant variations are ob- 
served) immediately after the magnetic field is introduced. 
To generate a one- dimensional pressure map, we fix the value 
of y so that a horizontal cut (in the x direction) along the 
upper surface of the computational domain passes through 
the strongest magnetic feature. The three plots in Fig. 3 
show the surface pressure distributions along such a cut at 
three different times. These plots clearly illustrate the sce- 
nario that was described in the previous paragraph. As the 
magnetic pressure grows, there is a corresponding decrease 
in the gas pressure as fluid rapidly drains out of the surface 
regions of the magnetic feature. Once the feature has formed 
(lower plot of Fig. 3), convective motions are strongly sup- 
pressed and the gas pressure in the magnetic region is very 
much smaller than that of the surrounding field- free fluid. 
Although partial evacuation, with accompanying field in- 
tensification, has already been observed in several previous 
studies (e.g. Hurlburt & Toomre 1988; Weiss et al. 1996; 
Vogler et al. 2005), the level of evacuation is much more 
dramatic in these simulations. 

The most surprising aspect of Fig. 3 is the fact that 
the magnetic pressure within the upper layers of the result- 
ing magnetic flux concentration is much larger than the gas 
pressure of the surrounding non-magnetic fluid. Since the 
ambient gas pressure increases rapidly with depth (P ga s oc 
(1 + 10z) 2 in the unperturbed polytropic atmosphere), the 
difference between the surrounding gas pressure and the in- 
ternal magnetic pressure decreases as we move away from 
the surface. For this magnetic feature, the magnetic pres- 
sure is comparable to the external gas pressure at a depth 
of approximately z = 0.05. However, even down to depths of 
approximately z = 0.1, the sum of the internal gas and mag- 
netic pressures is still larger than the external gas pressure. 
So, even below the surface, this magnetic feature is not sim- 
ply in pressure balance with its non-magnetic surroundings. 
Vogler et al. (2005) have shown that magnetic curvature 
effects can play a key role in the confinement of magnetic 
flux concentrations. However, in this case the field lines are 
predominantly vertical near the upper surface, so magnetic 
curvature effects are negligible. Note that the lack of field- 
line curvature also implies that these magnetic features are 
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Figure 3. The ^-dependence of the pressure distributions (each at fixed y) at three different time intervals along a horizontal cut through 
the strongest magnetic feature at the upper surface. Solid lines correspond to the magnetic pressure, P m ag; dashed lines correspond to 
the gas pressure, P ga s; the dotted lines represent the dynamic pressure, Pdyn ( see text). These snapshots are taken at (a) t=0.12, (b) 
t=0.61 and (c) t=1.61. 



not being confined in deeper layers (where the gas pressure 
is higher), since this would require the flux concentrations 
to spread out in the upper layers. 

In the absence of a confining influence due to magnetic 
curvature, the observed pressure imbalance implies that the 
surrounding convective motions must be playing an active 
role in the confinement of this magnetic feature. Put differ- 
ently, where the flow converges on an intense flux concentra- 
tion there has to be an excess stagnation pressure. Here this 
excess pressure is provided by the magnetic field in a region 
where the gas pressure is reduced. This is in total contrast 
with the simulation of umbral convection by Schiissler and 
Vogler (2006), where a divergent rising plume is associated 
with a weaker field, and the excess pressure results instead 
from an enhancement of density (associated with buoyancy 
braking). The significance of the dynamic pressure in these 
calculations is illustrated in Fig. 3 - it is clear that the dy- 
namic pressure around the magnetic feature can often be 
comparable to the local gas pressure. This shows that the 
dynamical influence of the surrounding convective motions 
cannot be ignored when considering models of photospheric 
magnetic field intensification. It should also be stressed that 
the magnetic energy density of the flux concentration that 



is illustrated in Fig. 3 is much larger than the mean kinetic 
energy density (i.e. Pd yn /2) of the surrounding granular con- 
vection. These simulations provide confirmation of the fact 
that the process of partial evacuation is important in the 
production of such super-equipartition fields. Without such 
a reduction in the local gas pressure, it is difficult to see 
how convective motions alone could produce such a strong 
magnetic feature. 

All the discussion so far has focused upon the early 
stages of these numerical simulations. The process of mag- 
netic field intensification is arguably the most interesting 
phase of the dynamics, although the evolution of the result- 
ing magnetic features also raises important issues. In the 
present idealised model, the net imposed vertical magnetic 
flux is a conserved quantity, so there will always be some 
accumulations of vertical magnetic flux somewhere within 
the computational domain. However, what is less clear is 
whether or not the occurrence of super-equipartition mag- 
netic fields is a transient feature of the model. Fig. 4 shows 
the magnetic field distribution several convective turnover 
times after the magnetic field was first introduced. It is 
clear that the most intense flux concentrations are still par- 
tially evacuated, and it is also clear that the magnetic en- 
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Figure 4. The magnetic field distribution at t = 26.26 (after sev- 
eral convective turnover times). Top: Like Fig. 2, this shows con- 
tours of the vertical magnetic field component superimposed upon 
the temperature profile in a horizontal plane at (z/d) = 0.05. Bot- 
tom: The pressure distributions along a horizontal cut through 
the strongest magnetic feature at the surface. As in Fig. 3, the 
solid line corresponds to the magnetic pressure, the dashed line 
corresponds to the gas pressure, whilst the dynamic pressure is 
represented by a dotted line. 

ergy density of these regions still exceeds the mean kinetic 
energy density of the surrounding convection. In addition, 
the fact that the internal magnetic pressure (in the upper 
layers) is still much larger than the external gas pressure 
indicates that the surrounding convective motions are still 
playing a key role in the confinement of the magnetic fea- 
ture. During their evolution, these features continuously in- 
teract with the convective motions, which implies that they 
are deformed, shredded and advected around the domain in 
a time-dependent fashion. However, despite these complex 
interactions, ultra-intense super-equipartition magnetic fea- 
tures seem to be a robust feature of the simulation and are 
not simply a transient phenomenon. 



3.3 The effects of varying £o 

One of the key parameters in these simulations is the mag- 
netic Reynolds number, Rm. Computational restrictions 
limit the range of values of Rm that can be considered, and 
all values that can be simulated will be very much smaller 
than real photospheric values. However, the simulation that 
has already been described qualitatively illustrates some of 
the main physical processes that occur during photospheric 
magnetic flux amplification. In this Section, we assess the 
effects of repeating this simulation for different values of £o- 
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Figure 5. The magnetic flux distribution during the kinematic 
phase for frj = 1.2. Like Fig. 2, this plot is for t = 0.24 and shows 
the spatial variation of the temperature and the vertical magnetic 
field component in a horizontal plane at (z/d) = 0.05 (i.e. near 
the upper surface). For ease of comparison, the same scales have 
been used in Fig. 2 and Fig. 5 for both the magnetic field and the 
temperature. 



Table 1. The (^-dependence of the magnitude of the vertical 
magnetic field component at a fixed position at the upper surface 
of the computational domain. The time is fixed at t = 0.12, and 
the centre of the magnetic feature is at (x, y) = (3.16,2.39) in 
all cases. Magnetic field strengths are normalised in terms of the 
strength of the imposed magnetic field. Also shown (in the middle 
column) is the corresponding value of Rm for each value of £o- 



Co 


Rm 


B z (max)/ B z (initial) 


0.2 


117 


3.61 


0.3 


83 


3.50 


0,1 


62 


3.40 


0.6 


42 


3.23 


1.2 


22 


2.85 


2.4 


11 


2.42 



This is equivalent to varying the magnetic Reynolds number 
of the flow, which is inversely proportional to (q: smaller val- 
ues of Rm correspond to larger values of £o and vice- versa. 
In order that comparisons can easily be made between the 
different cases, all magnetohydrodynamical simulations are 
started from exactly the same hydrodynamic initial condi- 
tions. 

Even during the brief kinematic phase, clear trends are 
observed as the parameter £o is varied. The effects of some 
of these trends are illustrated in Fig. [5] for C,o = 1-2. Like 
Fig. [2] this shows the spatial distribution of the temperature 
and the vertical component of the magnetic field near the 
upper surface of the computational domain. In order to allow 
direct comparison between Fig.[2]and Fig. [5] this snapshot of 
the simulation is also taken at t — 0.24, and the contours are 
plotted at the same levels - this implies that the scales for 
temperature and magnetic field correspond directly to those 
used in Fig. [2] In this kinematic phase, the temperature field 
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Figure 6. Probability density functions for the vertical component of the magnetic field at the upper surface of the computational 
domain. These functions correspond to different values of (g. These values are Co = 0.2 (top), Co = 0.6 (middle) and Co = 1.2 (bottom). 
For smaller values of Co (or, equivalcntly, larger values of Rm) there is a higher probability of stronger fields and there is a greater 
proportion of reversed polarity magnetic flux. 



is the same in both plots; however it is immediately apparent 
that the magnetic field structures in Fig. [5] are weaker and 
less localised than those seen in the Co = 0.2 case. This 
is a consequence of the increased importance of diffusion 
for larger values of Co (equivalently, lower values of Rm). 
The influence that Co has upon the magnetic field in this 
kinematic phase is quantified in Table [1] Here, we measure 
the strength of the vertical component of the magnetic field 
at (x,y,z) = (3.16,2.39,0) and t = 0.12, as a function of 
Co- This spatial location corresponds to the centre of the 
strongest surface magnetic feature at t — 0.12, which does 
not depend upon Co in the kinematic regime (although the 
location of this peak field is obviously a function of time). 
The table confirms that, at a fixed time, there are weaker 
fields at lower values of Rm. For simple flows, it is possible 
to construct a similarity solution for this kind of kinematic 
flux concentration, which leads to a power-law relationship 
between the peak magnetic field and the magnetic Reynolds 
number (see, e.g., Proctor & Weiss, 1982). Unfortunately, 
due to the complexity of the flow patterns in this case, there 
appears to be no analogous scaling here. 

For Co = 0.2, it was found that the high magnetic pres- 
sure in the strongest magnetic field concentrations quickly 



led to the partial evacuation of these magnetic features. 
Since these magnetic flux concentrations grow more gradu- 
ally for larger values of Co, this partial evacuation process is 
much slower in these cases and the effects of diffusion tend to 
play a more dominant role in limiting the flux concentration. 
Therefore, the amount of evacuation that occurs decreases 
with increasing Co- This demonstrates that the rapid inten- 
sification that occurred in the Co = 0.2 case is very much 
a high magnetic Reynolds number phenomenon. However, 
some evacuation is observed in all cases that were investi- 
gated except Co = 2.4, and wherever features do become 
partially evacuated, super-equipartition magnetic features 
are observed. Fig.|6]shows time-averaged probability density 
functions (pdfs) for the vertical component of the surface 
magnetic field, for three different values of Co- In all cases, 
the pdfs peak at B z = 0, which implies that the majority 
of the domain is field-free. In addition, there is a signifi- 
cant component of reversed magnetic flux at lower values of 
Co- At higher Rm there is a greater tendency for flux to be 
advected with any fluid flow, and these convective motions 
certainly have the potential to reverse magnetic flux at the 
granular boundaries at the surface of the domain. Most in- 
terestingly, there is not a large difference between the peak 
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fields in the pdfs for fo = 0.6 and £o = 0.2. Equivalently, 
the magnitude of the peak field appears to be only weakly- 
dependent upon Rm in this particular parameter regime. 
This suggests that rather than being controlled by diffusive 
effects, the peak attainable magnetic field in this parameter 
regime is determined by the combined effects of the exter- 
nal gas and dynamic pressures. So, although ultra-intense 
magnetic features form more rapidly at higher Rm, these 
features can also be produced at lower values of Rm, pro- 
vided that magnetic diffusion does not completely inhibit 
the flux concentration process before the magnetic field can 
become dynamically active. 



4 CONCLUSIONS 

This paper describes results from a series of numerical ex- 
periments that were designed to investigate the formation 
of localised magnetic flux concentrations at the solar photo- 
sphere. By adopting an idealised model, it was possible to 
assess the effects that varying the magnetic Reynolds num- 
ber might have upon this flux intensification process. As ex- 
pected, magnetic flux tends to accumulate preferentially in 
convective downflows, where it forms localised features - the 
horizontal scale of these features decreases with increasing 
magnetic Reynolds number. High magnetic pressures lead 
to the partial evacuation of these features as they form. At 
high values of Rm, the resulting field strengths are typically 
much larger than the equipartition value at which the local 
magnetic energy density balances the mean kinetic energy 
density of the surrounding granular convection. In addition, 
the strongest magnetic fields that are formed in the upper 
layers of the domain exert a magnetic pressure that is signif- 
icantly larger than the external gas pressure. The appear- 
ance of these ultra-intense magnetic fields shows that the 
dynamic pressure that is associated with the surrounding 
convection must be playing a key role in the confinement 
of these magnetic features to localised regions. Some super- 
equipartition magnetic fields are found in all cases except 
the most diffusive case (corresponding to £o = 2.4). 

It is interesting to relate these results to convective col- 
lapse models (see, e.g., Webb & Roberts 1978; Spruit & 
Zweibel 1979). Although those models are certainly an ide- 
alised representation of photospheric magnetic field inten- 
sification, our simulations suggest that they correctly iden- 
tify the most important process in the formation of super- 
cquipartition magnetic fields, namely the partial evacuation 
of magnetic regions by convective downflows. Our simula- 
tions do, however, raise some important issues relating to 
convective collapse. First, convective collapse models typi- 
cally adopt an initial condition that corresponds to a thin 
flux tube embedded in a static atmosphere. Since such mag- 
netic features form in well-established convective downflows, 
this static equilibrium never occurs in our simulations. Ad- 
ditionally, the evacuation process seems to begin before the 
flux accumulation process finishes. So our simulations could 
be seen as describing an "adjustment" rather than an insta- 
bility. The second point that is raised by these simulations 
concerns the pressure balance. Convective collapse models 
assume a constant balance between the external gas pressure 
and the gas and magnetic pressures within the flux tube. Our 
simulations indicate that the dynamical effects of the sur- 



rounding convection are playing a significant role in the con- 
finement of these magnetic features. The neglect of this dy- 
namic pressure underestimates the strength of the strongest 
magnetic fields that can be generated. Finally, Cameron & 
Galloway (2005) have modelled aspects of convective col- 
lapse by conducting numerical simulations of laminar mag- 
netoconvection in a simplified geometry. From these simula- 
tions they argue that super-equipartition magnetic features 
must be structured on a kinematic scale at the solar pho- 
tosphere. Although only a limited range of values for Rm 
has been considered here, our simulations do produce some 
super-equipartition fields on larger scales, and so contradict 
their view. 

When comparing results from these simulations with 
photospheric observations, it is important to remember some 
of the simplifying assumptions in this model. Processes such 
as radiative transfer have been neglected, so it is not possible 
to make detailed spectral comparisons between these simula- 
tions and observations. Nevertheless, this idealised approach 
has had a great deal of success in reproducing qualitative 
features of photospheric magnetoconvection. In qualitative 
terms, our simulations appear to be consistent with solar 
observations and provide a plausible explanation for the ap- 
pearance of super-equipartition magnetic features in the in- 
tergranular lanes in the quiet Sun. However, this model does 
have other limitations, notably the fact that many of the 
parameters that are used are not closely related to realistic 
solar values. We have also used highly idealised boundary 
conditions in these simulations: for example we impose a 
rigid boundary at the upper surface. It should be stressed 
that, although the transition to a subadiabatic stratification 
in the photosphere provides a "softer" boundary condition, 
this will not prevent the formation of extremely strong mag- 
netic fields. Another limitation that should be noted is the 
fact that (due to the periodic boundary conditions) the im- 
posed magnetic flux is independent of both z and t. This 
implies that net vertical magnetic flux cannot enter or leave 
the domain, therefore the initial non-zero magnetic flux is 
an invariant quantity in the model. In the quiet Sun, there 
is a continuous emergence of mixed polarity magnetic flux, 
which will interact with existing magnetic features. Such in- 
teractions will tend to limit the lifetimes of these magnetic 
features. This is something that cannot be represented in 
our idealised simulations, though it may be possible to make 
progress by careful choice of initial conditions. 

In the Introduction, we noted that it is still not clear 
whether quiet Sun magnetic features are simply fragments of 
reprocessed magnetic flux or whether they are generated lo- 
cally as the result of small-scale dynamo action. In fact, this 
model of compressible convection can drive a small-scale dy- 
namo once the magnetic Reynolds number exceeds a thresh- 
old of about Rm — 250. However, even for a marginally- 
excited dynamo, the magnetic Prandtl number is of order 2 
which is very much larger than the magnetic Prandtl num- 
ber in the solar photosphere. Whether or not such a dynamo 
could operate in the low magnetic Prandtl number regime 
has been the subject of considerable debate (see, for exam- 
ple, Boldyrev & Cattaneo 2004; Schekochihin et al. 2005). 
Although our idealised simulations do assume some pre- 
existing magnetic field (i.e. they do not generate this field 
self-consistently) , the processes of magnetic flux expulsion 
and intensification that are illustrated by these simulations 
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are generic and are therefore likely to be of relevance to the 
solar photosphere whether or not a local dynamo is operat- 
ing. 

This work is motivated by high resolution observations 
of the solar photosphere. Although current ground-based in- 
struments, such as the Swedish 1-metre Solar Telescope, are 
already provided very detailed images of the photosphere, it 
is likely that newer instruments (including those carried on 
the recently launched Hinode satellite) will reveal many new 
features of photospheric magnetoconvection. Over the next 
few years, these new observations will enable us to refine our 
current theoretical models, but will also inevitably present 
new theoretical challenges. 
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